Ordered vs disordered: Correlation lengths of 2D Potts models at /3 t 
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We performed Monte Carlo simulations of two-dimensional g-state Potts models with q = 10, 15, and 20 and 
measured the spin-spin correlation function at the first-order transition point fit in the disordered and ordered 
phase. Our results for the correlation length £d{fit) in the disordered phase are compatible with an analytic 
formula. Estimates of the correlation length £ {fit) in the ordered phase yield strong numerical evidence that 
R = Ufit)/Ufit) = l- 



1. INTRODUCTION 

Correlation lengths or inverse masses are 
among the most important quantities to charac- 
terize the properties of a statistical system. Un- 
fortunately there are only very few models for 
which the correlation length £ is exactly known 
and can thus serve as a testing ground for the 
employed numerical techniques. The best known 
example is the two-dimensional (2D) Ising model 
where £ is exactly known at all temperatures. But 
the situation is already much more involved for 
2D g-state Potts models with a partition function 
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l,...,g, (1) 



denote the lattice sites of a 



square lattice of size V = L x x L y , (ij) are 
nearest-neighbor pairs, and 5 SiS . is the Kronecker 
delta symbol. Here only the correlation length 
£d{Pt) in the disordered phase could be derived 
analytically [1], where fi t = \n(l + ^/q) is the first- 
order transition point of this model for q > 5. For 
the correlation length £ (/3*) in the ordered phase 
no analytical results are available. One goal of 
the present investigations was to test the conjec- 
ture that R = io{fit) I £d{fit) = 1/2, which was 
suggested quite heuristically by previous Monte 
Carlo (MC) studies [2]. 
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2. SIMULATION 

We studied the Potts model (1) in both the 
disordered and ordered phase and measured cor- 
relation functions at fi t for q = 10, 15 and 20 
on lattices of size V = L x L and V = 2L x L 
with L = 150,60 and 40 (« 14&). To take ad- 
vantage of translational invariance we used peri- 
odic boundary conditions. We carefully checked 
that our lattice sizes are large enough to suppress 
tunneling events such that, starting from a com- 
pletely random or ordered configuration, the sys- 
tem remained a sufficiently long time in the dis- 
ordered or ordered phase to perform statistically 
meaningful measurements. Analyses of autocor- 
relation times suggested to use in the disordered 
phase the single-cluster and in the ordered phase 
the heat-bath update algorithm. All error bars 
are estimated by means of the jack-knife tech- 
nique. 

To determine the correlation length £ we con- 



sidered the k 
tions, 
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of the spin-spin correlation function and its im- 
proved cluster estimator 



G(i,j) = (5 s 
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where Q(i,j) = 1, if i and j belong to the same 
Swendsen-Wang cluster, and 6 = otherwise. 
The projection removes the power-like prefactor 



in the large-distance behaviour of G and thus al- 
lows fits of j(")(i) = g {n) {i x , 0) to the Ansatz 

fl («)(x)= a ch(^/^) + 6ch(c^4^), (4) 
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where 



£<"> = C/^1 + (2tt</L 3/ ) 2 « C(l+ <n/L y ) 2 ). (5) 

In the ordered phase (4) can only be used for 
n^O. 



3. RESULTS 

3.1. Disordered phase 

A preliminary report of a first set of simulations 
in the disordered phase on L x L lattices was al- 
ready given in Ref.[3]. In the meantime we have 
further increased the statistics and added another 
set of simulations on asymmetric 2L x L lattices 
[4] . For both lattice geometries the first three en- 
ergy moments are found in very good agreement 
with the exact result for the average energy and 
with dual transformed large q expansions [4,5]. 

In the disordered phase we concentrated on the 
k y = projection of the correlation function (3). 
In a first step we fixed £(°) = £* at its theoretical 
value (see Table 1) and optimized only the re- 
maining three parameters in (4). For q = 10 the 
resulting fits to the L x L and 2L x L data are 
shown in Fig. 1 as dotted and solid lines, respec- 
tively. While the lines are excellent interpolations 
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Figure 1. Semi-log plot of the projected correla- 
tion functions for q = 10 in the disordered phase. 



Table 1 

Correlation length £d(Pt) in the disordered phase. 





q = 10 


q= 15 


g = 20 


Lx L 


9.0(5) 


3.70(16) 


2.24(6) 


(2L x L) x 


10.2(9) 


3.59(10) 


2.23(5) 


(2L x L) y 


9.3(7) 


3.62(16) 


2.33(7) 


exact 


10.5595... 


4.1809... 


2.6955... 



of the data over a wide range up to x « (5 . . . 6)£,d, 
we also see a clear tendency of the fits to lie sys- 
tematically above the data for very large x. In un- 
constrained fits to the Ansatz (4) this is reflected 
by a systematic trend to underestimate ^. By 
restricting the fit interval to larger x values the 
estimates for £* increase, but then also the sta- 
tistical errors grow rapidly [4] . Table 1 shows our 
results for £,d(Pt) using the fit range [x m - m ,L/2] 
with x m { n ss 2£d- We further estimate c« 1.5 — 2, 
independent of q, which stabilizes to c = 1.5(1) if 
£d is held fixed at its theoretical value. 

3.2. Ordered phase 

Also in the ordered phase we first checked that 
the average energy agrees with the exact result, 
and compared the higher moments with the large 
q expansions of Ref.[5]. For the specific heat see 
Table 2. Furthermore we looked at the magne- 
tization m = (q (m&x{rii})/V — l)/(q — 1) and 
its cluster estimator ml = (|C| max )/V, where 
rii denotes the number of spins of "orientation" 
i = 1, . . . , q and |C| max is the size of the largest 
(spanning) cluster. Also here we find very good 
agreement with the exact answer [6] (for q = 10 
and 2L x L, e.g., m = m! = 0.857 113(49), m ex = 
0.857106...). 

To determine £ we followed Gupta and Irback 
[2] and studied the k y = 2it / L y projection gW (x) 
of G. Since this removes constant background 



Table 2 

Specific heat at /3 t in the ordered phase. 



10 



15 



20 



Lx L 17.95(13) 
2Lx L 17.81(10) 
large q 18.1(1) 



8.016(21) 
8.004(19) 
8.00(3) 



5.351(15) 

5.3612(55) 

5.362(5) 
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Figure 2. Semi-log plot of the projected correla- 
tion function g^ for q = 10 in the ordered phase. 
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Figure 3. Ratio of effective correlation lengths in 
the ordered and disordered phase for q = 10. 



terms, g^(x) can be fitted with the Ansatz (4) 
also in the ordered phase. For q = 10 the qualita- 
tive behaviour of g^ is illustrated in the semi-log 
plot of Fig. 2. The comparison with g(°> of the 
disordered phase suggests that the two correla- 
tion functions are governed by the same asymp- 
totic decay law, i.e., that £ o (0t) = £d(Pt)- In fact, 
the dotted line interpolating the g^ data is a 
constrained fit to (4) assuming that £ = £<j = 

10.5595 To make this statement even more 

convincing we have plotted in Fig. 3 the ratio 
^ eff = tf/tf, where £ eff = l/\n[g(x)/g(x + 1)] 
is the usual effective correlation length (with the 
correction (5) for £„ already taken into account). 
The corresponding plots for q = 15 and q = 20 
look similar. 

4. DISCUSSION 

The numerical results for £d(Pt) confirm the an- 
alytical expression to about 10 — 20%. The sys- 
tematic underestimate is presumably caused by 
higher excitations which are neglected in the fits. 
Using a non-zero momentum projection in the or- 
dered phase we believe that £ (Pt) is of about 
the same accuracy. By comparing the two cor- 
relation length at (3 t we obtain strong numeri- 
cal evidence that £ = £,*. At first sight this is 
in striking disagreement with a very recent ex- 
act proof [7] of the earlier conjecture £ = £d/2 
for a certain definition of the ordered correlation 
length, £ 0j i. For another definition, £ 0j 2, how- 



ever, only the relation £ 0j2 > £ 0j i could be estab- 
lished, which would be consistent with our results 
if we identify the numerically determined £ with 
£o,2- We are currently investigating this problem 
in more detail [8] by using precisely the defini- 
tions of Ref.[7] which are based on geometrical 
properties of Potts model clusters such as, e.g., 
their diameter. 
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